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X-RAY SPECTROSCOPY OF THE ACCRETING MILLISECOND PULSAR XTE J0929-314 IN OUTBURST 
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ABSTRACT 

We present the high-resolution spectrum of the accretion-powered millisecond pulsar XTE J0929-314 during 
its 2002 outburst, measured using the Low Energy Transmission Grating Spectrometer onboard the Chandra X- 
ray Observatory . The 1.5-25.3 A (0.5-8.3 keV) Chandra spectrum is well fit by a power-law + blackbody model 
with photon index V = 1.55±0.03, blackbody temperature kT\,\, = 0.65±0.03 keV, and blackbody normalization 
^bb.km/^iokpc = 7.6 ±0.8. No emission or absorption features are found in the high-resolution spectrum, with 
a 3a equivalent width upper limit of < 0.007 A at 1.5 A and < 0.12 A at 24 A. The neutral absorption edge 
depths are consistent with the estimated interstellar absorption along the line of sight to the source. We found no 
orbital modulation of the 2-10 keV X-ray flux, to a 3er upper limit of 1.1%, which implies an upper limit on the 
binary inclination angle of i < 85° for a Roche-lobe-filling companion. We also present the broadband spectrum 
measured over the course of the outburst by the Rossi X-ray Timing Explorer (RXTE). The RXTE spectrum of 
XTE J0929-314 is also well fit with a power-law + blackbody model, with average values of T = 1.76±0.03, 
kTob = 0.66±0.06 keV, and /?bb,km/^iokpc = 5.9 ± 1.3 in the 2-50 keV energy range. The blackbody flux remained 
constant over the course of the outburst, while the power-law flux was strongly correlated to the (decreasing) flux 
of the source. We find that the difference in power-law photon indices measured from Chandra and RXTE spectra 
can be explained by a change in the power-law photon index at low energies. 

Subject headings: binaries: close — stars: neutron — pulsars: individual (XTE J0929-314) — X-rays: binaries 



1. INTRODUCTION 

Millisecond pulsars (MSPs) have long been considered one 
of the possible endpoints of low-mass X-ray binary (LMXB) 
evolution. The neutron star (NS) is thought to be spun-up to 
a millisecond period by accretion from its low-mass compan- 
ion. After the accretion phase has ended, the NS may turn on 
as a radio MSP. In the last five years, this theory has been con- 
firmed with the identification of three accretion-powered MSPs, 
SAX J1808.4-3658, XTE J1751-305 and XTE J0929-314 
(Wijnands & van der Klis 1998; Chakrabarty & Morgan 1998; 
Markwardt et al. 2002; Galloway et al. 2002). 

Interestingly, all three accretion-powered MSPs are in short 
period binaries, with the two most recently discovered having 
binary periods w43 min (Markwardt et al. 2002; Galloway et al. 
2002). These short periods place the MSPs XTE J175 1-305 
and XTE J0929-314 in the class of ultracompact binaries, de- 
fined as having orbital periods <80 min. Ultracompact bina- 
ries require hydrogen-deficient or degenerate donors (Nelson, 
Rappaport, & Joss 1986). Recently, O and Ne emission and ab- 
sorption features were discovered in two known and three sus- 
pected ultracompact systems (Schulz et al. 2001; Juett, Psaltis, 
& Chakrabarty 2001; Juett & Chakrabarty 2002). These results 
led the authors to conclude that the donor stars are degenerate 
C-O WDs. On the other hand, the XMM-Newton spectrum of 
XTE J1751-305 did not show emission or absorption features 
(Miller et al. 2002). In addition, no unusual abundances were 
required to fit the neutral absorption edges. This result is con- 
sistent with the suggestion of Markwardt et al. (2002) that the 
donor in XTE J175 1-305 is a He WD (see also, Bildsten 2002). 

XTE J0929-314 was discovered in April 2002 by the All 
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Sky Monitor onboard the Rossi X-ray Timing Explorer (RXTE; 
Remillard 2002). Further RXTE observations detected 185 Hz 
pulsations modulated by a 43.6-min binary orbit (Remillard, 
Swank, & Strohmayer 2002; Galloway et al. 2002). Radio and 
optical counterparts with positions consistent with the X-ray 
position were also detected (Rupen, Dhawan, & Mioduszewski 
2002; Greenhill, Giles, & Hill 2002; Cacella 2002). An opti- 
cal spectrum of XTE J0929-314 revealed emission lines from 
C III/N III A4640-4650 and Ha A6563 (Castro-Tirado et al. 
2002). Given the optical detection of emission lines, and the 
ultracompact nature and high-Galactic latitude of the source, 
XTE J0929-314 is an ideal candidate for a high-resolution X- 
ray spectroscopic study to search for emission and absorption 
features similar to those seen in other ultracompact systems. In 
this letter, we present results from a Director's Discretionary 
Time observation of XTE J0929-314 with the Chandra X-ray 
Observatory, as well as spectral results from the RXTE pointed 
observations throughout the outburst. 



2. observation and data reduction 

The 2002 outburst of XTE J0929-314 began around 
MJD 52370 and peaked at w 31 mCrab (2-10 keV) on 
MJD 52394 (Galloway et al. 2002). The pointed RXTE ob- 
servations commenced 2 d later, and proceeded once every few 
days through MJD 52456. The source dropped below the 3cr 
detection limit of 7.5 x 10~ 12 ergcm" 2 s _1 (2-10 keV; equiva- 
lent to 0.32 mCrab) after MJD 52443. 

XTE J0929-314 was observed throughout its 2002 outburst 
by the two pointed instruments onboard RXTE. The Propor- 
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tional Counter Array (PCA; Jahoda et al. 1996) consists of 
five gas-filled proportional counter units (PCUs) with a total 
effective area of w 6000 cm 2 , sensitive to X-ray photons in the 
2.5-60 keV range. The two instrument clusters comprising the 
High-Energy X-ray Timing Experiment (HEXTE; Gruber et al. 
1996) present an effective area of w 1600 cm 2 to photons in the 
energy range 15-250 keV. We extracted average spectra from 
standard instrument mode data ("Standard-2" and "Archive" 
for PCA and HEXTE, respectively) from each PCU/cluster for 
each observation. Since the instrumental gain is known to 
vary between PCUs, and also with time, we generated a sep- 
arate response matrix for each PCU and each observation us- 
ing pcarsp version 8.0, supplied with LHEASOFT version 
5.2. We estimated background spectra using the "combined" 
gain epoch 5 (beginning 2000 May 13) faint-source model as 
input to pcabackest version 3.0. Over the course of the out- 
burst, we collected 38 observations of XTE J0929-3 14 totaling 
123.5 ks. 

We also observed XTE J0929-314 with Chandra on 2002 
May 15 (MJD 52409) for 18 ks using the Low Energy Trans- 
mission Grating Spectrometer (LETGS) and the Advanced 
CCD Imaging Spectrometer (ACIS; Weisskopf et al. 2002). 
The LETGS spectra are imaged by ACIS, an array of six CCD 
detectors. The LETGS/ACIS combination provides both an 
undispersed (zeroth order) image and dispersed spectra from 
the grating with a first order wavelength range of 1.4-63 A 
(0.2-8.9 keV) and a spectral resolution of AA = 0.05 A. The 
various orders overlap and are sorted using the intrinsic energy 
resolution of the ACIS CCDs. The observation used a Y-offset 
of +1 .'5 in order to place the O-K absorption edge on the back- 
side illuminated S3 CCD, which has suffered less degradation 
than the front-side illuminated CCDs. 

Using tgdetect, we determined the zeroth order source 
position of XTE J0929-314 to be: R.A.=09 h 29 m 20?15 and 
Dec=-31°23'04."3, equinox J2000.0 (90% confidence error of 
0."6 3 ). The Chandra position is consistent with both the opti- 
cal and radio counterpart positions (Rupen et al. 2002; Green- 
hill et al. 2002; Cacella 2002). The first order dispersed 
spectrum of XTE J0929-314 had an average count rate of 
7.64±0.02 counts s" 1 . We examined the total count rate, as 
well as the count rates in two different energy ranges, to check 
for changes in the spectral state. We found no evidence for any 
change of spectral state during the Chandra observation. 

The "level 1" event file was processed using the CIAO 
v2.2 data analysis package 4 . The standard CIAO spectral 
reduction procedure was performed. We filtered the event 
file retaining those events tagged as afterglow events by 
the acis_detect_af terglow tool. Since order-sorting 
of grating spectra provides efficient rejection of background 
events, the afterglow detection tool is not necessary to detect 
cosmic ray afterglow events. No features were found that might 
be attributable to afterglow events. 

For bright sources, pileup can be a problem for CCD detec- 
tors (see, e.g., Davis 2002). The zeroth order Chandra spectra 
of XTE J0929-314 was heavily affected by pileup and was not 
used in this analysis. In addition, it was found that the first or- 
der spectrum of XTE J0929-314 suffered from pileup in the 
range 5-12 A (1-2.5 keV). In order to use the grating pileup 
kernel in ISIS, we created detector response files (ARFs) for 
each chip using the CIAO tool mkgarf and then combined 

3 See http://asc.harvard.edu/cal/ASPECT/celmon/index.html 
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them to create +1 and -1 order ARFs using a custom tool de- 
veloped by J. Davis of the HETGS instrument team. (This tool 
is similar to the standard CIAO tool but correctly calculates the 
fractional exposure at each response bin. This functionality has 
been added to the standard tools in CIAO v2.3.) The grating 
pileup kernel models the effect of pileup on grating spectra, 
similar to the pileup model available for CCD spectra in ISIS 
and XSPEC (see, Davis 2002, for a discussion of grating pileup 
modeling). In addition, the data and responses were rebinned 
to 0.083 A to reflect the size of an ACIS event detection cell 
(3 CCD pixels). Background spectra were extracted from the 
standard LETG background regions. Spectral analysis of the 
Chandra observation of XTE J0929-314 was performed using 
ISIS (Houck & DeNicola 2000). 



3. SEARCH FOR ORBITAL MODULATION 

The presence (or absence) of orbital modulation of the X-ray 
flux can help constrain the inclination of the binary, and in turn 
the mass of the companion. Galloway et al. (2002) determined 
the orbital period of XTE J0929-314 to be 43.6 min from the 
Doppler modulation of the pulse arrival times. We searched 
both the RXTE and Chandra data for orbital modulation of the 
X-ray flux. 

For the RXTE analysis, we calculated background- 
subtracted, solar-system barycenter corrected, 16-s binned 
lightcurves for each observation from photons in the energy 
range 3-13 keV. This range was selected to maximize the 
signal-to-noise ratio, and is identical to that used for the pulse 
timing analysis of Galloway et al. (2002). For each observa- 
tion, we created a phase folded lightcurve of 8 bins using the 
measured orbital parameters, with phase zero set to the epoch 
of 90° mean longitude, T n n. The length of each observation 
was short compared to the decay timescale of the X-ray flux of 
XTE J0929-314, so it was not necessary to correct for the de- 
cline in X-ray flux during a single observation. Each lightcurve 
was then fit with a sine curve with variable amplitude and phase. 
From these results, we place a 3cr upper limit on the fractional 
rms modulation of 1 . 1 %. 

We also searched the Chandra data for orbital modulation of 
the X-ray flux. For the Chandra analysis, we first barycentered 
and randomized the first order grating events. Randomizing 
of the event arrival times consists of adding a random quan- 
tity uniformly distributed between 0-3.2 s, in order to avoid 
aliasing caused by the readout time. We calculated the phase 
of each first order event using the orbital parameters from Gal- 
loway et al. (2002) and created an 8 bin phase-folded lightcurve. 
We found an upper limit of 2.3% rms, consistent with the RXTE 
results. 

Since the Chandra data is not subject to time gaps due to 
spacecraft motion, we also calculated a Fourier transform of 
the first order events which had been made into a 10-s binned 
lightcurve. We searched for modulations of the X-ray flux with 
frequencies between 6xl0" 5 and 5xl0" 2 Hz. We found no 
evidence for any periodic modulation, with a 90%-confidence 
upper limit of 1 .4% for the fractional rms amplitude over the 
frequency range. The 90%-confidence upper limit for the frac- 
tional rms amplitude of a signal at the orbital period is 0.4%, 
well within the upper limits given by the phase folding analy- 
ses. 
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FIG. 1 . — Spectral fit parameters derived from fitting an absorbed blackbody 
+ power-law model to RXTE observations of XTE J0929-314. The panels, 
from top to bottom, show the blackbody temperature r bb in keV; the black- 
body normalization (= Skm/^lOkpc). estimated bolometric flux of the black- 
body component (1(T 10 ergcirT 2 s -1 ), power-law photon index T, and power- 
law normalization (photons keV -1 cnT 2 s -1 at 1 keV). The solid line shows the 
weighted mean over all the observations, while the dashed line shows the ±la 
limits. The dotted line shows the time of the Chandra observation. Error bars 
show the Ict uncertainties. 



4. SPECTRAL ANALYSIS 

4. 1 . RXTE Spectral Analysis 

We fitted the combined 2-25 keV PCA and 15-50 keV 
HEXTE spectra with a number of models which typically give 
good fits for other X-ray pulsars, including broken and cut- 
off power laws and analytical Comptonization approximations. 
There was no evidence for a high-energy spectral cutoff within 
the energy range in which the source was detected (typically 
< 50 keV). We did however measure significant residuals be- 
low 10 keV. These residuals could be minimized by adopting a 
broken power law model where the spectral index decreased by 
f=s 10% (i.e. the spectrum became harder) above 7 keV, or by 
adding a blackbody component with kT\$, w 0.5 keV to a power 
law. While both these models gave fits of similar quality, we 
rejected the former as unphysical and instead adopted the latter 
for the remainder of the RXTE fits. 

We also included a multiplicative component to take into ac- 
count the attenuation by intervening (neutral) material with cos- 
mic abundances. The equivalent hydrogen column density Nh 
exhibited no significant variations between observations, and in 
the mean was <~ 10 21 cm" 2 , consistent with the line-of-sight 
values interpolated from dust and H I maps (7.6 x 10 20 and 
10 21 cm -2 respectively; Schlegel, Finkbeiner, & Davis 1998; 
Dickey & Lockman 1990). For consistency with the fits to 
Chandra data, we froze the Nh value in subsequent fits at 



7.6 x 10 cm . The mean value of the reduced-% (= xt) f° r 

5 For more information see, http://cxc.harvard.edu/cal/Links/Acis/acis/Cal_prods/qeDeg/index.html 



power-law + blackbody fits to spectra from all the observations 
was 1.04. While the worst-fitting observation gave a % 2 = 2.44, 
the largest residuals arose from variations between PCUs (par- 
ticularly at the lower energy bound) rather than any broad trend 
with energy. By excluding these lowest channels from the fit, 
and assuming a systematic error of 1%, we were able to reduce 
the fit statistic for that observation to xt =1-30. 

The weighted mean of the blackbody temperature, kT^, 
was 0.66 ± 0.06 keV, while the normalization v^i* = 
^bh.km/^iokpc = 5.9 ±1.3 (all errors are quoted at 90%- 
confidence unless otherwise noted). Both parameters showed 
variations during the outburst; however, they were significantly 
anticorrelated, so that the estimated bolometric blackbody flux 
(oc XbbTbJ,) remained constant to within our measurement un- 
certainties (x 2 = 1.06, see Figure 1). In the first observation, 
the blackbody component comprises only 3% of the flux in the 
2-60 keV band, whereas in the last observation where signifi- 
cant blackbody and power-law components were detected, the 
blackbody accounted for 16% of the flux. 

We calculated a weighted mean photon index of T = 1 .76 ± 
0.03. The principal source of flux variation was the normal- 
ization of the power-law component (see Figure 1). While 
all parameters exhibited variations throughout the outburst, 
the power-law normalization was the only parameter to show 
a significant correlation with the integrated 2-10 keV flux 
(Spearman's rank correlation coefficient p = 0.908, equiv- 
alent to 5 Act; Press et al. 1992). A power-law compo- 
nent was not significantly detected in the last observation on 
MJD 52443. A blackbody-only fit for that day exhibits a signifi- 
cantly (6.6cr) lower integrated blackbody flux than the weighted 
mean level for the previous observations. We did not detect 
XTE J0929-3 14 in subsequent observations (3a upper limit on 
the 2-10 keV flux of 7.5 x 10" 12 erg cm" 2 s" 1 ). 

4.2. Chandra Spectral Analysis 

We fit the Chandra spectrum with a power-law + blackbody 
model including absorption. The +1.'5 Y-offset in the obser- 
vational setup allows for the entire 1.5-25.3 A (0.5-8.3 keV) 
+1 order spectrum to be imaged on the S3 CCD. This CCD is 
backside-illuminated and has suffered less degradation than the 
front-side illuminated CCDs. Because of this and the lack of 
counts above 25 A, we performed our spectral fits on only the 
+1 order spectrum in the range 1.5-25.3 A. 

Of the three known X-ray MSPs, XTE J0929-314 has the 
highest Galactic latitude and, probably the smallest interstel- 
lar column density along the line of sight to the source. Our 
Chandra observation allows for a direct determination of the 
strength of the absorption edges, in particular the O-K edge and 
the Fe-L triplet. To account for absorption and measure the 
optical depths of the edges, we used the tbvarabs absorp- 
tion model (see, Wilms, Allen, & McCray 2000) with O and Fe 
abundances set to zero and the A^h fixed to 7.6 x 10 20 cm" 2 , the 
hydrogen column from dust maps (Schlegel et al. 1998). The O 
edge was fit with an edge model, while the Fe-L edges were 
fit using a custom model that employs the optical constant mea- 
surements of Kortright & Kim (2000). In addition, the C abun- 
dance in tbvarabs was allowed to vary in order to account 
for the known instrumental contamination of the ACIS CCDs 5 
(Plucinsky et al. 2002). We also included Gaussian lines to fit 
the interstellar atomic O absorption line at 23.5 A and the ab- 
sorption feature 



4 



JUETT, GALLOWAY, & CHAKRABARTY 




2 4 6 8 10 12 

Wavelength (A) 




Wovelength (A) 

FIG. 2. — (upper panels) Chandra LETG +1 order flux spectrum (in units 
of photons cm" 2 s" 1 A" 1 ) of XTE J0929-3 14 with best-fit power-law + black- 
body model with absorption. The arrows indicate the positions of the O edge, 
and the atomic O and Fe2C>3 absorption lines, (lower panels) Fractional resid- 
uals ([data-model]/model) of the LETG spectral fit shown above. The absorp- 
tion models are consistent with the predicted interstellar absorption from dust 
maps. 

at 23.36 A, attributed to a ls-2p transition in Fe2C>3 (see, Schulz 
et al. 2002). 

We found a best-fit photon index of 1.55±0.03, with a nor- 
malization of (4.04±0.15)xl0" 2 photons keV" 1 cm" 2 s" 1 at 1 
keV. The Chandra best-fit blackbody temperature was kT\fo = 
0.65 ± 0.03 keV, with # b b,kmMok P c = 7.6 ± 0.8. Our best-fit 
model had a xt = 1001 and an unabsorbed 2-10 keV flux of 
2.7xl0" 1() ergcm" 2 s" 1 . 

The best-fit optical depth for O was 0.34±0.07, which 
translates to an equivalent hydrogen column density of 
(1.2±0.3)x 10 21 cm" 2 , using the cross-section of Henke, Gul- 
likson, & Davis (1993) and the O ISM abundance of Wilms 
et al. (2000). Although this is somewhat higher than the mea- 
sured hydrogen column from dust maps, an instrumental con- 
tribution to the O edge has been reported, with an optical depth 
of w0.10 (Plucinsky et al. 2002). This contribution would add 
Nu ~ 0.36 x 10 21 cm" 2 to the expected value, making the Chan- 
dra O edge measurement compatible with the dust map mea- 
surement. The best-fit line wavelength and equivalent width 
(EW) of the atomic O absorption line were 23.52±0.05 A and 
0.053^o q3 5 A. The Fe203 absorption line and Fe-L edges were 
not significantly detected, with an upper limit on the Fe column 
density of 3 x 10 16 cm" 2 , equivalent to anNn < 1.1 x 10 21 cm" 2 . 
The best-fit C abundance was 7.8+q 2 times the interstellar C/H 
ratio of Wilms et al. (2000). This is slightly lower than the ex- 
pected instrumental + interstellar contribution of 9.9±0.5 times 
the interstellar ratio. The instrumental edge depths for C and O 
were estimated using the recent calibration results presented at 



the Chandra X-ray Center website 5 (see also, Plucinsky et al. 
2002). 

Besides the interstellar O absorption line, there were no 
prominent emission or absorption lines in the high-resolution 
spectrum of XTE J0929-314 (see Figure 2). We performed a 
careful search of the Chandra spectral residuals to place lim- 
its on the presence of any spectral features. Gaussian models 
with fixed FWHM = 2000 km s" 1 were fit at each point in the 
wavelength range 1 .5-25.3 A. From this, we can place a 3cr up- 
per limit on the EW of any line feature, either emission or ab- 
sorption. The EW limit increases approximately linearly with 
wavelength, varying from 0.007 A at 1.5 A to 0.12 A at 24 A. 

Our measured C abundance depends on the assumed contin- 
uum model since the edge depth can not be directly measured 
(the C edge at 43 A is outside our bandpass). In order to test if 
the C abundance was the cause of the difference in the photon 
index between the Chandra and RXTE fits, we refit the Chandra 
spectrum with a fixed C abundance of 9.9 times the interstellar 
C/H ratio. With the C abundance fixed, the best-fit photon in- 
dex was 1.62±0.03. This is still significantly different from the 
best-fit photon index found in the RXTE spectra. The best-fit 
power-law normalization increased to (4.50±0.07)x 10~ 2 pho- 
tons keV" 1 cm" 2 s" 1 . The other fit parameters were consistent 
with the previous fit within errors. We also considered the pos- 
sibility that the pileup model was giving rise to the lower photon 
index in the Chandra spectral fits. When the Chandra data is fit 
without the pileup model, and excluding the piled-up region (1- 
3 keV) of the spectrum, we find spectral parameters consistent 
with the results found using the pileup model. 



5. COMBINED Chandra AND RXTE SPECTRAL FITS 

To investigate the difference in the spectral results be- 
tween Chandra and RXTE, we performed combined spectral 
fits with the Chandra and RXTE data. A single RXTE ob- 
servation (70096-03-05-00), taken within a few hours of the 
Chandra data, was used to represent the RXTE spectrum of 
XTE J0929-314. To simplify fitting, the piled-up region (1- 
3 keV) of the Chandra data was excluded and the pileup ker- 
nel was not used. The combined data were fit with a black- 
body + either a power-law, broken power-law or the comptt 
model which describes Comptonization of soft photons by a hot 
plasma (Titarchuk 1994). The broken power-law and Comp- 
tonization models allow for turnover at low energies which 
could explain the difference in the Chandra and RXTE spectral 
results. The absorption was described by the tbvarabs model 
with the C abundance allowed to vary for the Chandra observa- 
tion only, while best-fit edge depths from the Chandra analysis 
were used for the O and Fe edges. We also included a Gaussian 
line, fixed to the best-fit Chandra value, to model the atomic O 
absorption at 23.5 A. Multiplicative constants were included to 
compensate for any instrumental normalization differences. We 
find that the normalizations of the PCA spectra are «20% larger 
than for Chandra, while the HEXTE spectra are «20% lower. 
Similar normalization differences have been noted before (see, 
e.g., Kuulkers et al. 2002, and references therein). 

As a check on the instrumental calibrations, we directly com- 
pared the Chandra and RXTE data in the same energy range by 
fitting the high energy (3-8 keV) Chandra data, in conjunction 
with the RXTE data, to an absorbed power-law + blackbody 
model with absorption fixed to 7Y H = 7.6 x 10 20 cm" 2 . The fits 
are consistent with the results of the RXTE data, T = 1 .74 ± 
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0.02, kT bh = 0.66 ± 0.03 keV, and /? b b,kmMokpc = 5.9 ± 0.7, 
with xj = 1.03. 

We then fit the full Chandra and RXTE spectra with a power- 
law + blackbody model. When the edge parameters were fixed 
to the Chandra best-fit values, the photon index had a best- 
fit value of 1.628±0.009 and with a xl = 1-20 (see Table 1). 
We note that the xl of the power-law + blackbody fit to the 
RXTE data alone was 1.06. If we fix the C abundance to the ex- 
pected value of 9.9 times the interstellar C/H ratio, the photon 
index increases to 1.656±0.009 (xl = 1.14). If we instead al- 
low the C abundance to vary, the best-fit photon index increases 
to 1 .72±0.02 (xl = 1 -09), giving a photon index consistent with 
the RXTE spectral results, but with a C abundance of 14.1 ±1.2 
times the interstellar C/H ratio. This value is inconsistent with 
the estimate of the C abundance from instrument calibration 
work (see §4.2). As can be seen, the continuum spectral pa- 
rameters are highly dependent on the assumed absorption. 

We performed the same fits using the comptt in place of 
the power-law component. The Comptonization model pro- 
vides for a low-energy turnover in the spectrum. The electron 
temperature was fixed to 50 keV, the upper limit of the energy 
range. In addition, we used a spherical geometry in order to 
compare with the spectral fit of SAX J1808.4-3658 (Titarchuk, 
Cui, Wood 2002). The comptt model provides a better fit to 
the data compared with the power-law (xl = 1.09-1.07, see Ta- 
ble 1). We performed an F-test to formally compare the two 
models and found that in all cases the Comptonization model 
was better with a significance > 99%. For the two fixed values 
of the C abundance, the results of the comptt + blackbody 
fits are consistent within error with a seed photon temperature, 
kTo w 0.13 keV and optical depth, t p w 2.7. When the C abun- 
dance is allowed to vary, the fit improves but becomes inconsis- 
tent with the other results, which is probably due to degeneracy 
between the parameters. 

In the best-fit comptt + blackbody model, the black- 
body component has parameter values of kT bb = 0.71 keV and 
^bb.km/^iokpc = 5.5. These values are consistent with the power- 
law fits both to the combined data and with the blackbody pa- 
rameters from each instrument fit, indicating that the black- 
body temperature and flux are independent of the exact con- 
tinuum model used. We find kTo < kT bb , similar to the fit of 



SAX J1808.4-3658 using the same model (Titarchuk et al. 
2002). We note that the comptt model is in itself an ap- 
proximation to a full treatment of Comptonization, in partic- 
ular at soft energies. If we assume that the blackbody provides 
the input photons to the Comptonization component, then we 
would expect kTo = kT bb . We fit the data with the constraint 
kT = kT bb , but the fit was unreasonable with xl > 2.0. In- 
terestingly, Gierlihski, Done, & Barret (2002) used a differ- 
ent Comptonization model and found kTo = kT bb m spectral fits 
of SAX J1808.4-3658. We employed the same model to fit 
XTE J0929-314 (compPS; Poutanen & Svensson 1996), but 
found similar values for the input photon temperature and opti- 
cal depth as found using the comptt model. 

We also fit the data with a broken power-law + blackbody 
model to test for low-energy turnover without relying on the 
assumptions inherent in the Comptonization model. We find 
that the data is also well fit (xl = 108) by a broken power- 
law + blackbody model with Ti = 1.26-1.56, T 2 = 1.75, and 
Fbreak = 1 -4—4A keV, and blackbody parameters consistent with 
previous fits (see Table 1). 

6. DISCUSSION 

We have found that the spectrum of XTE J0929-314 is 
well described by a power-law + blackbody model over lim- 
ited energy ranges, with interstellar absorption consistent with 
Nh = 7.6 x 10 20 cm" 2 , the expected hydrogen column density 
along the line of sight to the source. The Chandra spectrum 
shows no prominent emission or absorption lines in the spec- 
trum (similar to XTE J175 1-305; Miller et al. 2002), with an 
EW limit which varies from < 0.007 A at 1.5 A to < 0.12 A 
at 24 A. No orbital modulation of the X-ray flux was detected 
in either the Chandra or RXTE data. From the lack of X-ray 
eclipses, we set an upper limit on the binary inclination of 
i < 85° for a Roche-lobe-filling companion. In addition, the 
absence of dips in the X-ray lightcurve suggests that the source 
is most likely at an inclination of i < 60° -70° (see, e.g., Frank, 
King, & Lasota 1987; White, Nagase, & Parmar 1995). 

The blackbody spectral component in LMXBs is generally 
attributed to the neutron star, possibly from a "hot spot" where 
the accretion column meets the neutron star surface. From our 



Table 1 



Combined Chandra and RXTE Spectral Fits 



Model 


Carbon 
Abundance b 


(keV) 


£break 

(keV) 


r 2 /r 


A^po/A'comp 


(keV) 


-ftkm/dlOkpc 


X 2 /dof 


PL 
PL 
PL 


7.8 (fix) 

9.9 (fix) 
14.1±1.2 


1.628±0.009 
1.656±0.009 
1.72±0.02 






4.58±0.07 
4.90±0.08 
5.7±0.3 


0.730±0.015 
0.717±0.016 
0.68±0.02 


5.3±0.3 
5.3±0.3 
5.6±0.4 


1006/837 
950/837 
904/833 


Comp 
Comp 
Comp 


7.8 (fix) 

9.9 (fix) 
13.2±1.2 


0.144±0.009 
0.125±0.014 

0.07!?- 04 




2.72±0.06 
2.69±0.06 
2.55±0.07 


3.21±0.09 
3.6±0.2 

5.83:§ 


0.709±0.016 
0.707±0.017 
0.67±0.02 


5.4±0.3 
5.5±0.4 
5.8±0.4 


911/836 
906/836 
894/832 


BKNPL 
BKNPL 


7.8 (fix) 

9.9 (fix) 


1.26±0.09 
1.56±0.03 


1.4±0.3 
4.4±0.9 


1.74±0.02 
1.75±0.02 


5.06±0.19 
4.67±0.17 


0.66±0.03 
0.61±0.03 


5.8±0.7 
7.7±1.9 


900/835 
897/835 



"All errors are quoted at the 90%-confidence level 
b Carbon abundance ratio relative to ISM, (C/H)/(C/H) K m 

c Power-law and broken power-law normalization at 1 keV in units of 10~ 2 photons keV -1 cnf 2 s -1 . Comptonization 
normalization in units of 10~ 3 . 
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fits, we find a blackbody radius of /?bb.km = (5.3-7.7) c/iok pc - 
Given the lower limit on the distance of 5 kpc (Galloway 
et al. 2002), we find a lower limit on the blackbody radius of 
Rhb > 2.7 km. We note that this lower limit does not include 
a correction for the conversion between the color temperature 
and the effective temperature of the neutron star atmosphere 
(see, Lewin, van Paradijs, & Taam 1993, for a discussion). This 
correction could increase the inferred radius by a factor of m 2. 
But even after this correction, the inferred radius is not consis- 
tent with canonical NS models for a distance of 5 kpc suggest- 
ing a "hot spot" emission region. 

During the outburst, we found that the blackbody flux re- 
mained constant, while the power-law flux declined. In addi- 
tion, when the combined Chandra and RXTE spectrum is fit 
with a Comptonization model, we find that the blackbody tem- 
perature, fcTbb, is significantly greater than the seed photon tem- 
perature of the Comptonization model, kTo. This evidence sug- 
gests that the power-law and blackbody components are associ- 
ated with different emission processes. 

We compared our results with two physical models that were 
developed to explain the spectral properties of the accretion- 
powered MSP SAX J1808.4-3658 (Gierlihski et al. 2002; 
Titarchuk et al. 2002). Both models were formulated using data 
from the 1998 outburst of SAX J1808.4-3658, but come to 
very different conclusions regarding the relationship between 
the power-law or Comptonization component and the black- 
body. Gierlihski et al. (2002) assumed that the blackbody was 
the source of the seed photons (i.e, kTo = kT\,\,) and proposed that 
the Comptonization component comes from the shock heated 
accretion column with soft photon input from the "hot spot" 
blackbody component. Such a model implies that the black- 
body and power-law/Comptonization component fluxes should 
be correlated, which was the case for SAX J1808.4-3658 but 
was not seen for XTE J0929-3 14. 

On the other hand, Titarchuk et al. (2002) found that kT < < 
A:7t,b and suggested that the accretion disk provided the major- 
ity of the soft photon input to a spherical Compton cloud, of 
radius R w 60 km, from which the Comptonization component 
originates. The Titarchuk et al. (2002) model allows for the the 
blackbody and power-law/Comptonization component to vary 
independently. Titarchuk et al. (2002) also suggested that the 
low (to ~ 4) optical depth of the Comptonizing region found 
for SAX J1808.4-3658 allows for the pulsations of the NS to 
be detected, while the higher optical depths (to > 5) reported 
in more luminous LMXBs suppress the pulsation amplitudes. 
We find a best-fit t w 2.7 for XTE J0929-3 14, consistent with 
this explanation. However, we note that Miller et al. (2002) did 
not find an acceptable fit of the XMM spectra of the millisecond 
X-ray pulsar XTE J175 1—305 using a Comptonization model. 
While this model seems more appropriate for the spectral re- 
sults of XTE J0929-314, Titarchuk et al. (2002) do not discuss 
what emission component is responsible for the pulsations. 

It is important to realize that both models are based on data 
with a low-energy bound of 2 keV, which is above the peak flux 
of the thermal components. The lower energy range of Chan- 



dra allows us to make a more robust measurement of kTo for 
XTE J0929-314 than was possible for the RXTE spectrum of 
SAX J1808.4-3658. We conclude that neither model is com- 
pletely appropriate to explain both our results. More obser- 
vations of these sources throughout outburst are necessary in 
order to fully understand their phase-averaged spectral proper- 
ties. We suggest that in-depth data analysis, including pulse- 
phase spectroscopy, will provide more insight into the relation- 
ship between the various spectral components. The pulse-phase 
resolved spectral analysis of XTE J0929-3 14 will be presented 
in a later paper. 

We find that the best-fit photon index was significantly dif- 
ferent between the Chandra and RXTE fits with T = 1.55-1 .62 
for the Chandra fits, dependent on the assumed instrumental 
C abundance, and T = 1.76 for the RXTE fits. Since fits with 
and without the pileup model give consistent parameter values, 
we reject the possibility that pileup causes the lower Chandra 
photon index. The combined Chandra and RXTE spectrum of 
XTE J0929-3 14 is well fit by either a Comptonization + black- 
body model (which implicitly includes a low-energy turnover 
from the Comptonization component), or a broken power-law 
+ blackbody model, with break energy 1 .\-\A keV. The high- 
energy (E > 3 keV) Chandra spectrum is consistent with the 
RXTE spectral results, indicating that instrumental differences 
alone do not give rise to the spectral turnover. Based on these 
results, we suggest that the difference in best-fit photon indices 
between Chandra and RXTE arises from spectral turnover at 
low energies. While a power-law model is a reasonable ap- 
proximation to Comptonization at high energies, extrapolating 
the power law to low energies leads to a predicted flux that di- 
verges at zero energy. Therefore, the power-law must turnover 
at some energy for the integrated to flux to be finite. 

In addition, the difference in the photon index between the 
Chandra and RXTE spectral fits is similar to that noted by 
Miller et al. (2002) between the RXTE and XMM spectral fits of 
XTE J 1 75 1 -305 . From their XMM EPIC spectrum, Miller et al. 
(2002) found a best-fit power-law photon index of 1 .44, while 
Markwardt et al. (2002) found a best-fit photon index of 1 .7-1 .9 
using RXTE. We interpret the difference between the XMM and 
RXTE spectral results for XTE J175 1-305 as evidence of low 
energy turnover in the broadband spectrum of that source, pro- 
viding independent support that the turnover is astrophysical in 
origin. 

We would like to thank Chandra Director Harvey Tanan- 
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ful discussions with Craig Markwardt, Jon Miller, Mike 
Nowak, and Dimitrios Psaltis. This research has made use 
of data obtained through the High Energy Astrophysics Sci- 
ence Archive Research Center Online Service, provided by the 
NASA/Goddard Space Flight Center. This work was supported 
in part by NASA under contract NAS8-01 129 and grant NAG5- 
9184. 
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